Remarks on Cosmic String Formation during Preheating on Lattice Simulations 
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We reconsider the formation of (global) cosmic strings during and after preheating by calculating 
the dynamics of a scalar field on both two- and three-dimensional lattices. We have found that 
there is little differences between the results in two and three dimensions about the dynamics of 
fluctuations, at least, during preheating. Practically, it is difficult to determine whether long cosmic 
strings which may affect the later evolution of the universe could ever be produced from the results of 
simulations on three-dimensional lattices with smaller box sizes than the horizon. Therefore, using 
two-dimensional lattices with large box size, we have found that cosmic strings with the breaking 
scale r, ~ 10 16 GeV are produced for broad range of parameter space in rj, while for higher breaking 
scales (?j~3x 10 16 GeV), their production depends crucially on the value of the breaking scale rj in 
our simulations. 
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One of the hot topics at the reheating stage after infla- 
tion is the possibility of the formation of topological de- 
fects . This phenomenon is due to large nonthermal 
fluctuations during preheating (6[j^] and efficient rescat- 
tering 0,|ll[ both of which are caused by the Bose en- 
hancement effects. 

At the preheating stage, very large nonthermal fluc- 
tuations are produced, (<5</> 2 ) ~ c 2 M 2 where M p is the 
Planck mass and c = 1CT 2 - 1CT 3 ||g]. These fluctu- 
ations make the shape of the effective potential of the 
field <fi change to the type of the potential which has a 
minimum at the origin if the potential V((f>) is of spon- 
taneous symmetry-breaking type. It may be regarded 
that the symmetry is restored Later, when the 

amplitude of these fluctuations is redshifted away by the 
cosmic expansion, the symmetry is spontaneously bro- 
ken and topological defects may be created. Thus, the 
mechanism for producing the topological defects seems 
to be somewhat similar to the Kibble mechanism in high 
temperature theory. 

For the order estimation of the critical value of the 
breaking scale where cosmic strings are not formed 
above that value, it is sufficient to see the amplitude of 
the nonthermal fluctuations produced during preheating: 
(64> 2 ) 1/2 ~ 10 16 GeV 0. This order estimation is in good 
agreement with numerical estimations done in Refs. [Q,^). 
But numerical simulations on the lattices which follow 
the dynamics of the scalar fields reveals several new re- 
sults, such as the effects of the rescattering or the 
number of defects estimation B . 

In our previous paper B, we investigated the dynam- 
ics of a complex scalar field using two-dimensional lattice 
simulations, taking into account that the box size is large 
enough to cover the horizon size and that the lattice size 
is small enough to identify cosmic strings safely, and con- 
cluded that (long) cosmic strings would not be produced 



if the breaking scale r\ was larger than r\ ~ 3 x 10 16 GeV. 
On the other hand, the authors of Ref. [|5j showed the 
possibility that cosmic strings could be formed even if 
rj = 6 X 10 16 GeV using the three-dimensional lattice 
where the box size is smaller than the horizon. In this 
paper, we will discuss both two- and three-dimensional 
lattice simulations, and make certain that these results 
are consistent, commenting on the limitations of both 
simulations. 

First we show that results from lattice simulations in 
two and three dimensions are not different from each 
other when we study the parametric resonance during 
preheating. To be concrete, let us consider a complex 
scalar field with the effective potential: 
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where A is a small coupling constant. This model has 
a global U(l) symmetry, and cosmic strings are formed 
when the symmetry is spontaneously broken. What we 
have to do is to integrate the equation of motion: 
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For numerical simulations, it is convenient to use rescaled 
variables: 
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where $ = |^(0)|- Setting a(0) = 1, we obtain 
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where fj = i]/$o and the prime denotes differentiation 
with respect to r. The second term of the left-handed 
side can be omitted since the energy density of the uni- 
verse behaves like radiation at the early times and also 
the scale factor becomes very large later. In this case, 
the rescaled Hubble parameter, h (r) = H(t)/\/\&o, and 
the scale factor a(r) become [Q 

h(r) = ^a- 2 (r), (7) 

and 

«W = yT + l, (8) 

respectively, when $ is assumed to be an inflaton (even 
if $ is not an inflaton, results are the same as in the case 
of rescaling the breaking scale in an appropriate way, see 
Ref. Q) and we set a(0) = 1. For the initial conditions 
we take 

x ee Re^(O) = 1 + 5x(x), (9) 
y ee Im^(O) = 5y(x), (10) 

where the homogeneous part comes from its definition 
(we call x direction for real direction and y for imagi- 
nary), and 5x, y(x) is a small random variable of 0(1O -7 ) 
representing the fluctuations. We also set small random 
values for velocities. 

Since the physical length and the horizon grow pro- 
portional to a and a 2 , respectively, the rescaled horizon 
grows proportional to a. The initial length of the hori- 
zon is 4(0) = 3/a/2 « 2.12. Therefore, the rescaled 
horizon size grows as £h( T ) = 4(0)a(r). It is thus better 
to take the box size larger than £(r en d) where T enc i is the 
time at the end of the simulation. This is because those 
strings whose lengths are shorter than the horizon scale 
will be in the form of loops, which will shrink and dis- 
appear very soon. Only longer strings than the horizon 
will survive to affect the later evolution of the universe. 
On the other hand, the width of the topological defect is 
(VAf?) -1 which corresponds to (^a(r)) -1 in the rescaled 
units. Since it decreases with time, one lattice length 
should be at least comparable with the defect width at 
the end of the calculation. 

Leaving these facts aside, let us just compare the 
evolution of fluctuations on two-dimensional lattices 
with three-dimensional ones. Here we take 128 3 three- 
dimensional lattices and 4096 2 two-dimensional lattices 
with the lattice size A£ = 0.3 for both cases. We find no 
difference between both growth exponents (see Fig. [l]). 
Notice that the data used on the two-dimensional lat- 
tice (top panel) is the same used in Fig. 4 of Ref. j|, 
where it was plotted linearly instead of logarithmically 
as in Fig. 111. f\ Moreover, we obtain very similar results 



The authors of Ref. M might be misleaded by this point, 
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FIG. 1. Evolution of the fluctuations for 77 = 3 x 10 16 GeV 
on two- (top panel) and three- (bottom panel) dimensional 
lattices. The solid and dotted line denote the fluctuations in 
y and x directions, respectively. Notice that the data used 
on the two-dimensional lattices is the same used in Fig.4 of 
Ref. [4]. 

in two- and three-dimensional lattices, at least, on the 
effects of parametric resonance during preheating. 

Since small loop strings will shrink and disappear very 
soon, and cause no influence on the universe, we are inter- 
ested only in infinitely long strings. In three-dimensional 
lattices, they can be considered as those strings that pen- 
etrate through the box of the lattice (come into the box 
from one side and go out to the other side) and survive 
until the later time. Actually, we have found tempo- 



since they claimed that the growth exponent of the fluctuation 
in the x direction is much larger in three-dimensional lattices 
than that in the two-dimensional cases of our previous results 
in Ref. §. 
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rary formation of cosmic strings at any breaking scale, 
which confirms the result of Ref. [0. In most values 
of the breaking scale, however, these strings are in the 
form of small loops, and disappear very quickly. We 
have also found strings longer than the box size. In 
Fig. B, we take the breaking scale 77 = 3.08 x 10 16 GeV, 
the lattice size A£ = 0.3, and integrate the equation of 
motion until r = 280. If we integrated until the time 
before r = 265, we would conclude that cosmic strings 
were formed for this breaking scale. However, as we can 
see in Fig. [| these long strings feel the attractive force 
from each other, and they form into loops beyond the 
box size, which is much smaller than the horizon scale: 
NA£ = 38.4 < 4(r = 280) » 282, where N is the 
number of the lattice. 

We agree with the authors of Ref. || that longer strings 
tend to be created when the velocity of the real part 
of the oscillating homogeneous mode x becomes almost 
zero at the moment when it passes through the origin 
of the effective potential (in their words, when the mo- 
ment of the symmetry breaking nearly coincides with the 
moment when (</>i)(= x) passes through zero |5j). They 
argued that it leads to the fact that the long string for- 
mation is a non-monotonic function of the breaking scale 
0, which we confirm to some extent. However, in more 
than one hundred runs of our simulations, for almost all 
the breaking scales higher than 3 x 10 16 GeV, we have 
not found long strings stretched out beyond the box size 
except for two breaking scales: 77 ~ 3.08 x 10 16 GeV and 
77 ~ 3.16 x 10 16 GeV on the 128 3 lattices with A£ = 0.3 
(Actually, long strings cannot be found at these scale, 
if another initial configurations of initial fluctuations are 
used). Even in these scales, long strings deform into loops 
and shrink and disappear very soon. 

Therefore, we cannot make any definite conclusion on 
the formation of cosmic strings in the sense that we can- 
not tell whether or not they may affect the evolution of 
the universe. This is because the box size cannot be taken 
to be larger than the horizon scale in three-dimensional 
lattices. We have thus studied the formation of cosmic 
strings in two-dimensional lattices in the previous paper 
P|, abandoning one dimension in space because of the 
lack of memory capacity of computers. 



*We have recently found long strings formed at 3.01 x 
10 16 GeV, but they are very unstable in the sense that strings 
in the form of loops develop into the form of long strings, and 
later they again deform into loops and disappear. Notice that 
the range of the breaking scale where long strings are formed 
is also very narrow at this scale as in the above two scales. Ac- 
tually, these scales correspond to the circumstances that the 
velocity of the real part of the oscillating homogeneous mode 
x becomes almost (but not exactly) zero at the moment when 
it passes through the origin of the effective potential, so there 
may be more breaking scales where long cosmic strings are 
formed. 







FIG. 2. Evolution of cosmic strings for 77 = 3.08 x 10 16 
GeV. We plot from r = 264 (top) to r = 270 (bottom). 
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FIG. 3. Dependence of the lifetime (jd — Tf) when long 
strings penetrate the box of the lattice on the breaking scale 77 
around 77 = 3.078 x 10 16 GeV. Here Td is the destruction time 
when a long string deforms into a loop, and 77 ~ 250 is the 
formation time when a long string is first formed. Notice that 
only loop strings can be formed outside the narrow strip of the 
breaking scale (77 < 3.075 x 10 16 GeV, 3.082 x 10 16 GeV < rf). 
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FIG. 4. Dependence of the lifetime (t<j — 77) when long 
strings penetrate the box of the lattice on the breaking scale 
77 around rj = 3.158 x 10 16 GeV. Notice that the formation 
time is 77 ~ 250, and only loop strings can be formed outside 
the narrow strip of the breaking scale (77 < 3.156 x 10 16 GeV, 
3.164 x 10 16 GeV < rj). 



Nevertheless, we can extract useful information from 
the results on three-dimensional lattices. If we take the 
criterion of the formation of (infinitely) long strings as 
the cases when those strings stretched across the box size 
have been once formed, for the sake of the conservative 
discussion, we have found that the corresponding cases 
are in the only very narrow region of the parameter space 
of the breaking scale for the fixed initial conditions: | 
A77/77 ~ 3x 1CT 3 . See Figs. || and f|. Figures || and || show 
the lifetime (r<j — 77) with respect to the breaking scale r) 
for - 3.08 x 10 16 GeV and - 3.16 x 10 16 GeV, respectively. 
We can see that the size of the narrow regions are both 
A 7/ ~ 0.009 x 10 16 GeV. Here r d is the time when the 
long strings destruct into loops, and 77 is the formation 
time for those strings. 

However, this criterion may not be so sufficient, since 
those strings stretched across the lattice box will deform 
into loops and disappear very soon. We can take other 
criterion, which may reflect the idea that it is important 
to consider a long cosmic string. When we observe the 
lifetime (r^ — 77), we find that long strings survive a few 
times longer at a certain range of the breaking scale, as 
shown in Figs. || and 0. We can expect that long strings 



^For other initial conditions, narrow ranges of the breaking 
scale where long cosmic strings are produced appear in the 
different scales, and they are also very narrow. 



that stretch beyond the horizon size will be formed only 
within such very narrow regions. In this case, therefore, 
the breaking scale should be in the very narrow ranges 
(A77/77 ~ 10 ) in order for (infinitely) long strings to 
be formed for the fixed initial conditions in our simu- 
lations (These features can be also seen on the lattices 
with larger box size: N = 200). In other words, the long 
string formation is very sensitive to the breaking scale. 

These results can be understood as follows. If there is 
no gradient force, the dynamics of the field is determined 
only by the homogeneous mode. But the initial values 
of the field at each point on the lattice differs from each 
other by O(10 -7 ), we naively expect that cosmic strings 
are formed only in the very narrow ranges of the breaking 
scale of O(10~ 7 ). Owing to the preheating stage, fluctu- 
ations become large, so that the degree of the narrow 
ranges is somewhat relaxed. However, it seems that the 
full development of rescattering does not occur yet for 
r\ ~ 3 x 10 16 GeV, as we mentioned in Ref. j|. We will 
see later that the fundamentally different feature appears 
for the lower breaking scales. 

As mentioned above, it is difficult to tell whether long 
cosmic strings are formed on three-dimensional lattices 
with a smaller box size than the horizon volume, and 
this is why we have calculated in two dimensions, in 
order to make definite conclusions using both two and 
three-dimensional simulations complementarily. Similar 
results are found in two-dimensional lattices. In the two- 
dimensional case, we cannot distinguish long strings from 
loops, because all the strings are assumed to be infinitely 
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FIG. 5. Dependence of the number of strings per hori- 
zon on the breaking scale r\ around r\ — 3 x 10 16 GeV at 
t = 410 (dashed), r = 500 (dotted), and r = 600 (solid), 
in two-dimensional simulations. 

stretched along the z direction. Instead, we observe the 
number of cosmic strings in the horizon size at each time. 
The physically meaningful criterion for the formation of 
long cosmic strings is whether the number of strings per 
horizon remain (almost) constant or not as time goes on. 
Even if we cannot distinguish very long strings from loop 
ones, we can regard very nearly located string-antistring 
pair in two dimensions as a small loop string in three di- 
mensions. Actually, they annihilate very soon, similar to 
small loop strings which will shrink and disappear very 
quickly. Thus, when we observe the number of strings 
as time goes on, it will remain at a certain value if long 
strings are formed, since loop strings (string-antistring 
pairs in terms of two dimensions) will disappear and very 
long strings (isolated strings in terms of two dimensions) 
will remain. Notice that we include both strings and an- 
tistrings in the numbers. As a result, we find that it 
depends on the breaking scale very crucially at the scale 
r\ ~ 3 x 10 16 GeV, as is seen in Fig. [|. At the scales 
X] ~ 3.02 x 10 16 GeV and r) ~ 3.09 x 10 16 GeV, a dozen 
cosmic strings are formed, and the numbers do not de- 
crease so much. On the other hand, at the other scales, 
the numbers of strings decreases, and finally, we find no 
strings at all. Thus, the value of the breaking scale has 
to be in the very narrow ranges such as the degree of 
Arj/r/ ~ 10~ 2 , where A77 is defined as the breaking scales 
at which the number of strings per horizon remains at 
a certain value as time goes on. This implies, together 
with the results in three dimensions, that long string for- 
mation is very sensitive to the breaking scale at these 
scales. 

Moreover, we have simulated 30 realizations of ini- 
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FIG. 6. Average numbers of strings per horizon at 
r = 600 over 30 realizations for each breaking scale rj around 
t) — 3 x 10 16 GeV in two-dimensional simulations. 

tial conditions for fluctuations for each breaking scale, 
and find that the dependence of the average numbers 
of strings per horizon on the breaking scales, shown in 
Fig. ||, coincides with the above particular one of Fig. ||. 
This not only confirms the crucial sensitiveness to the 
breaking scale, but also implies that the main factor 
which determines whether long strings are produced is 
the value of the breaking scale. Initial fluctuations do 
not affect the dynamics of the scalar field so much at 
i)~3x 10 16 GeV @. 

On the contrary to the cases with somewhat higher 
breaking scale such as 77 ~ 3 x 10 16 GeV, many cosmic 
strings are formed at 77 ~ 10 16 GeV as is seen in Fig. [?]. 
We have found more than a dozen strings per horizon size 
at any value of the breaking scales near 77 ~ 10 16 GeV, 
and the numbers do not decrease so much. We thus con- 
clude that the formation of cosmic strings occurs in the 
very broad region of the breaking scale at 77 ~ 10 16 GeV. 

Finally, we should comment on the relation between 
the symmetry restoration and the topological defect (cos- 
mic string) formation. To know whether the symmetry is 
restored or not is a difficult task to do, and its methods 
are somewhat uncertain. The authors of Ref. [gj argued 
it in terms of the shape of the effective potential for the 
scalar field. If the potential has a minimum at the origin, 
the symmetry is restored. We may make sure that the 
potential has a minimum at the origin in the following 
way, as done in Ref. ||. The field <j> (radial direction of 
$) will oscillate around the origin in the potential of the 
form (<p 2 — rj 2 ) 2 when its amplitude is larger than yplr\. 
But, as we can see in Fig. ||, cj) is oscillating around the 
origin even when its amplitude becomes smaller than the 
breaking scale T). It is thus very clear that the effective 
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FIG. 7. Dependence of the number of strings per horizon 
on the breaking scale r\ around 77 = 10 16 GeV at r = 1420 (dot- 
ted), and t — 1800 (solid) in two-dimensional simulations. 

potential should have a minimum at the origin. We can 
see this phenomena in all cases in Fig. ^. However, in the 
above, we actually see that topological defects are pro- 
duced only when the breaking scale is 77 = 10 16 GeV, not 
in the cases of 77 = 3 x 10 16 GeV and 77 = 6 x 10 16 GeV. 
Therefore, the symmetry must be fully restored only in 
the case of 77 = 10 16 GeV, where rescatterings play a cru- 
cial role for that. 

Then what is the relation between the symmetry 
restoration and defect production? We conclude that de- 
fect formation is the signal of the full restoration of sym- 
metry. We thus discriminate between symmetry restora- 
tion and the shape of the effective potential which has a 
minimum at the origin. In other words, you cannot tell 
if the symmetry is restored only by observing the shape 
of the potential. 

In conclusion, we have reconsidered the formation of 
(global) cosmic strings during and after preheating by 
calculating the dynamics of the scalar field on both two- 
and three-dimensional lattices, and confirmed the results 
of both Ref . || and our previous ones [Q . We have found 
that there are little differences between the results in two 
and three dimensions, at least, at the preheating stage. 
It is obvious that phase-space volume is larger in the 
three dimensions than two, and this effect might affect 
somehow in the rescattering stage, but we expect that 
it seems subdominant when we obverse our numerical 
simulations. 

Practically, it is difficult to determine whether long 
cosmic strings which may affect the later evolution of 
the universe could ever be produced from the results of 
simulations on three-dimensional lattices, since they will 
deform into large loops and disappear very soon because 
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FIG. 8. Evolution of the $-field for various breaking 
scales. 77 = 6 x 10 16 GeV, 77 = 3 x 10 16 GeV, and 77 = 10 16 GeV 
from the top to the bottom. Here the amplitudes are plotted 
in the physical unit, not in the rescaled one. 
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of the small box size of the lattices. Moreover, we have 
found that cosmic strings with a higher breaking scale 
than 3 x 10 16 GeV could only be produced in the very 
narrow ranges of the breaking scale in our simulations. In 
Rcf . M , they referred to this formation of cosmic string as 
the formation of a nonmonotonic function of the breaking 
scale. We confirm this result and, in addition, find that 
the formation of long cosmic strings occurs in very small 
parameter space of the breaking scale 77. In other words, 
it is very sensitive to the value of the breaking scale. In 
two-dimensional simulations, long strings and loops can 
be distinguished to some extent, even though it is diffi- 
cult, by observing the evolution of the number of strings 
per horizon, and we have found similar phenomenon that 
long strings are produced only within very narrow range 
of the breaking scale around 77 ~ 3 x 10 16 GeV for fixed 
initial conditions. On the contrary, they are produced for 
a wide range of the breaking scale when 77 ~ 10 16 GeV. 

Cosmologically, we are interested only in (infinitely) 
long cosmic strings stretched beyond the horizon. Those 
strings with breaking scale 77 < 10 16 GeV are naturally 
formed after preheating, since they are produced inde- 
pendent to the actual values of the breaking scale. In 
other words, they are produced and their numbers re- 
main constant in every horizon volume. 

On the other hand, when the breaking scale is larger 
(77 ~ 3 x 10 16 GeV), it is very difficult to connect our 
results directly to the actual probability of the string for- 
mation. What we have found is that the formation of a 
cosmic string with 77 ~ 3 x 10 16 GeV depends crucially 
on the breaking scale. As mentioned, since initial fluctu- 
ations do not affect whether long strings are formed so 
much, we can tell how many long strings are produced, if 
the value of the breaking scale and the initial condition 
for the homogeneous mode are fixed, which is, in general, 
determined if the inflation model is specified. 
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